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TECHNICAL MEMORANDUM 


IRRADIANCE FROM STARS BASED ON BLACKBODY THEORY 

INTRODUCTION 


This report describes a computer program which calculates a theoretical blackbody curve for a 
star when given only the star’s temperature and visual magnitude. The program was written to produce 
approximate stellar spectra to be used in the design of the Space Telescope Fine Guidance System (FGS). 
In addition to the blackbody curve, the program calculates the number of photons that would be received 
from the star by the Space Telescope per 25 msec interval in wavelength increments of 0.2 microns. 

The equations used for the computation work best for stars that emit primarily in the range from near 
infrared to ultraviolet. Cooler stars may not give accurate results. 


DESCRIPTION OF PROGRAM CALCULATIONS 


This computer program is designed to give data for a blackbody curve and a count of photons 
emitted by the star. The only inputs are the temperature and visual magnitude of the star being used. 

The first step in the program is to find the visible irradiance for a star [I(m y )] using equation (1). 


m v = -2.5 log jo [I(m v )/I Q ] 


(D 1 


Where 

I Q = irradiance of zero magnitude star 
m y - visible magnitude 

Equation (1) easily converts to equation (2) with I(m y ) dependent on I Q and m y . 
I(m y ) = I Q x 10 ^ mv ^ Watts cm~^ 

To find the blackbody visible response fraction [? 7 e (T)] equation (3) is used. 


r? e (T)=- 


/W X (T) Se x dX 


( 2 ) 


(3) 1 


JW X (T) dX 



Where 


Se = the human eye response^ 

W X (T) = the Planck function, equation (4) 

T = temperature in degrees Kelvin 
X = wavelength 

and the Planck function wavelength (X) is given by, 

(4) 3 

where 

c = speed of light 
k = Planck’s constant 
h = Boltzman’s constant 

The top integral in equation (3) is calculated in the program over lengths of 0.34 to 0.82 with a 
dX (delta wavelength) of 0.02 microns. The bottom integral is the Stefan Boltzman function which is 
equal to 

5.679 x 10“^ x T^ watts, cm^ (5)1 

To find the maximum value of the Planck function [W x max(T)J equation (6) is used. 

W x max(T) = 1.25 x 10"^ x T^‘ watts/cm^/M (6)1 

Using equations (2), (3), (5) and (6), the peak of the blackbody curve becomes equation (7). 

I(m v ) W x Max(T) j 

H-vpeak = x watts/cm^/u (7) 

A %(T) JW X (T) dX 


W X (T) = 


2ir c^ h 

X.5 (ehcAkt-i) 


Watts/cm^/p 


2 



From the value for H^peak and using equation (8), 


H(X,T) 


0.2900 

(XT) 5 (e 1 438 A T _i) 


x H^peak watts/cnA/p 


(8) 1 


values for the blackbody curve are calculated for each (wavelength) used. This data then can be used 
to create a plot for the blackbody star. 

To find the amount of watts/ju received, the blackbody data was multiplied by a collecting aper- 
ture (A). The aperture used in the sample is from the Space Telescope. To give an approximation in 
watts over an area around the wavelength increments the watts/ju is multiplied by 0.02 microns. 


H s j(X,T) = H(X,T) x A x 0.02 watts 


(9) 


Watts are divided by hc/X to give photons. As is shown in equation (10). 


H st (X,T) 

Hp(X,T) = — - — - — photons/sec 
hc/X 


( 10 ) 


CONCLUSION 


By the use of this computer program, theoretical values are created for a blackbody curve. These 
values are calculated from a star’s visual magnitude and may be used to produce photon counts at certain 
wavelengths for use in many applications. These applications include throughput analysis of lenses and 
mirrors where efficiency is important. 

The blackbody curve created with this program is only applicable outside the atmosphere. For 
operation inside the atmosphere, values of spectral absorption should be applied. The values calculated 
are assumed to be from a star which behaves as a “Planckian Emitter.” There is evidence that while many 
stars behave this way, some do not. This is especially true for cooler stars. 
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CONSTANTS 


*l Q = 3.1 x 10 - * 3 watts cm 2 
c = 2.99793 x 10^° cm/sec 
h = 6.626176 x 10~ 34 J.S 
k = 1.380662 x 10~ 23 J/K° 
n = 3.1415926 

A = 3.89 x 10 4 cm 2 for Space Telescope 


* Allen’s Astro physical Journal states that I Q = 3.6 x 10 * 3 watts cm 2 , but to be conservative 
the Applied Optics value of 3.1 x 10~* 3 watts cm 2 is used. 



EXAMPLE RUN OF PROGRAM 


This program is currently being used in a throughput analysis program for the Fine Guidance 
System of the Space Telescope. The values for H(\,T) are normalized (Fig. 1) and multiplied with the 
throughput value of the other surfaces (Fig. 2), producing a sample of the amount of photons that will 
reach the optical detectors that produce a signal for fine pointing of the Space Telescope. Two example 
runs of the program show the number of photons emitted by the star. 


MAGNETUDE= 13.50 TEMPERATURE= 5500/10 



Figure 1. Normalized blackbody curve of a star. 
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S P SYSTEM CONDITION ■ 

1.00 1 OTA TRANSMISSION - IRD SPEC 

1.00 II VH? MIRROR 

1.00 20 VLR SURFACE 

1.00 1 CLEAR FILTER 

1.00 1 TRANSMISSION OF ALL GLASS AM3 HOMOSIL 

1.00 1 BEAM SPLITTER 

1.00 1 QUANTUM EFFICIENCY OF D E TECTOR 

1.00 1 EFFICIENCY OF ELECTRONICS 



WAVELENGTH (MICRONS) 


Figure 2. OTA-FGS throughput curve. 
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SAMPLE RUN OF PROGRAM BEDATAPRQG 


CONTROL COMMANDS 
.'SET F : 20/EYERESP 
! SET F: TT ME; DRC? BIN 

START OF PROGRAM 
! START BBDATAPROG 


IN THIii BLACK BODY DATA PROGRAM, YOU WILL 
ENTER A MAGNETUDE AND TEMPERATURE <K) 

AND THE PROGRAM WILL COMPUTE THE PHOTONS PER 
WAVELENGTH. 


ENTER THE MAGNE'flJDI: AND TEMPERATURE <K> 
? 1 3 . 5 , 5500 



WAVELENGTH® 

.. 340 

H( LAMBDA)® 

.6855 HE- 17 

PHOTONS® 

223. 

WAVELENGTH® 

.. 360 

H< LAMBDA)® 

.7S983SE-17 

PHOTONS® 

278. 

WAVELENGTH™ 

.. 380 

H< LAMBDA)® 

. 833652E— 17 

PHOTONS® 

329. 

WAVELENGTH® 

..400 

H< LAMBDA)® 

. 964913E-17 

PHOTONS® 

378. 

WAVELENGTH-" 

.420 

H< LAMBDA)® 

. 103264E-16 

PHOTONS® 

425. 

WAVELENGTH® 

.. 440 

H( LAMBDA)® 

. 10S668E-16 

PHOTONS® 

468. 

WAVELENGTH® 

.460 

H < LAMBDA ) = 

. 1.12750E— 16 

PHOTONS® 

508. 

WAVELENGTH® 

.480 

H< LAMBDA)® 

. 115597E-16 

PHOTONS® 

543. 

WAVELENGTH® 

.500 

H<LAMBDA)= 

. 1 17323E-16 

PHOTONS® 

574. 

WAVELENGTH® 

. 520 

H( LAMBDA)® 

. 1 LS055E-16 

PHOTONS® 

601. 

WAVELENGTH® 

.540 

H( LAMBDA)® 

. 117921E-16 

PHOTONS® 

623. 

WAVELENGTH® 

.. 560 

H (LAMBDA)® 

. 117050E-16 

PHOTONS® 

642. 

WAVELENGTH® 

. 580 

H< LAMBDA)® 

. 1 1556* IE— 16 

PHOTONS® 

656. 

WAVELENGTH® 

.600 

H( LAMBDA)® 

. U3563E-16 

PHOTONS® 

667. 

WAVELENGTH® 

. 620 

HC LAMBDA)® 

. 111156E-16 

PHOTONS® 

675. 

WAVELENGTH® 

. 640 

H< LAMBDA)® 

. 108428E-16 

PHOTONS® 

679. 

WAVELENGTH® 

. 660 

H < LAMBDA ) = 

. 105454E— 16 

PHOTONS® 

681 . 

WAVELENGTH® 

. 680 

H( LAMBDA)® 

. 102302E— 16 

PHOTONS® 

631. 

WAVELENGTH® 

..700 

HI LAMBDA)® 

. 990262E— 17 

PHOTONS® 

679. 

WAVELENGTH® 

.. 720 

H< LAMBDA)® 

. 956753E-1 7 

PHOTONS® 

674. 

WAVELENGTH® 

..740 

H< LAMBDA)® 

.. 922884E- 1 7 

PHOTONS® 

669. 

WAVELENGTH® 

.. 7<!.0 

H< LAMBDA)® 

.68898 IE- 17 

PHOTONS® 

662. 

WAVELENGTH® 

. 780 

H< LAMBDA)® 

. 855313E— 17 

PHOTONS® 

653. 

WAVELENGTH® 

. 800 

H( LAMBDA)® 

. 822094E— 17 

PHOTONS® 

644. 

WAVELENGTH® 

. 820 

Ht LAMBDA)® 

.. 789495E— 1 7 

PHOTONS® 

634. 


*STOP* END OF PROGRAM 
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SAMPLE RUN OF PROGRAM BBDATAPROG 
Run #2 


CONTROL COMMANDS 


.'SET F s 20/EYERESP 
! SET Fs-TT ME? DRC? BIN 


START OF PROGRAM 


! START BBDATAPROG 


IN THIS BLACKBODY DATA PROGRAM, YOU WILL 
ENTER A MAGNETLIDE AND TEMPERATURE ( K ) 

AND THE PROGRAM HILL COMPUTE THE PHOTONS PER 
WAVELENGTH. 


ENTER THE MAGNETUDE AND TEMPERATURE (K> 
?14. 5. 5500 



WAVELENGTH- 

. 340 

H( LAMBDA) 55 

WAVELENGTH 55 

.. 360 

H( LAMBDA) 55 

WAVELENGTH- 

• J30 

H< LAMBDA) 55 

WAVELENCTH= 

.. 400 

H( LAMBDA) 55 

WAVELENGTH 5 " 

. 4:20 

H( LAMBDA) 55 

WAVELENGTH 55 

.440 

H (LAMBDA) = 

WAVELENGTH 55 

. 460 

H ( LAMBDA ) 55 

WAVELENGTH 55 

. 430 

H( LAMBDA) 55 

WAVELENGTH 55 

. 500 

H( LAMBDA) 55 

WAVELENGTH- 

. 520 

H( LAMBDA) 55 

WAVELENGTH™ 

.540 

H< LAMBDA) 55 

WAVELENGTH= 

.. 560 

H( LAMBDA) 55 

WAVELENGTH 55 

.580 

H< LAMBDA) 55 

WAVELENGTH® 5 

.600 

H( LAMBDA) 55 

WAVELENGTH 55 : 

. 620 

H( LAMBDA) 55 

WAVELENGTH- 

.640 

H( LAMBDA) 55 

WAVELENGTH 55 

. 660 

H( LAMBDA) 55 

WAVELENGTH 55 

.680 

H( LAMBDA) 55 

WAVELENGTH 55 

.. 700 

H< LAMBDA) 55 

WAVELENGTH 55 

. 720 

H< LAMBDA) 55 

WAVELENGTH 55 

.. 740 

H< LAMBDA) 55 

WAVELENGTH 55 

. 760 

H ( LAMBDA ) 555 

WAVELENEiTH 55 

. 780 

H< LAMBDA) 55 

WAVELENGTH 55 

..300 

H( LAMBDA) 55 

WAVELENGTH 55 

. 820 

H( LAMBDA) 55 

*STOP* END OF 

PROGRAM 


272903E-17 

PHOTONS 55 

91. 

31 444 IE -17 

PHOTONS 55 

111. 

351789E— 17 

PHOTONS 55 

131 . 

384140E— 17 

PHOTONS 55 

150. 

41 1 103E--17 

PHOTONS 55 

169. 

43261 7E- 17 

PHOTONS 55 

1 86. 

448867E-17 

PHOTONS 55 

202. 

46020 IE- 17 

PHOTONS 55 

216. 

467073E-17 

PHOTONS 55 

229. 

469986E-17 

PHOTONS 55 

239 . 

469455E— 17 

PHOTONS 55 

243. 

465937E--17 

PHOTONS 55 

256. 

460057E— 17 

PHOTONS 55 

261 . 

452105E— 17 

PHOTONS 55 

266. 

442523E-17 

PHOTONS 55 

269. 

431660E-17 

PHOTONS 55 

270. 

4 1 9823E- 1 7 

PHOTONS 55 

271. 

407272E—17 

PHOTONS 55 

271. 

394232E-1 7 

PHOTONS 55 

270. 

3S0892E— 17 

PHOTONS 55 

269. 

367408E-17 

PHOTONS 55 

266 . 

35391 IE- 17 

PHOTONS 55 

263. 

340507E— 17 

PHOTONS 55 

260. 

327283E-17 

PHOTONS 55 

256. 

314305E— 17 

PHOTONS— 

252. 
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*TYl-42 

1.000 C 

2.000 C 

3.000 C 

4.000 C 

5.000 C 

6.000 C 

7.000 C 

8.000 C 

9.000 
10. 000 

1 1 . 000 

12. 000 

13. 000 

14.000 

15. 000 

16.000 

17. 000 

18.000 

19.000 C 

20.000 C 

21.000 C 

22.000 

23. 000 

24.000 

25.000 

26.000 C 

27.000 C 

28.000 0 

29.000 

30 . 000 

31.000 

32.000 

33.000 

34. 000 

35.000 C 

36.000 C 

37.000 C 

33.000 

39.000 C 

40.000 C 

41.000 C 

42.000 

* 


*■*••»•*«•# BLACKBODY DATA PROGRAM 


THIS PROGRAM IS DESIGNED TO CREATE DATA FOR 
A BLACK BODY CURVE FROM A GIVEN TEMPERATURE 
AND A GIVEN MAGNETUDE 
WRITTEN BY WILLIAM JACOBS 

DIMENSION S< 100) , W ( 100 ) , R ( 100 ) , PHOTON < 100 ) 

D I MENS I ON HST < .1 00 ) , H ( 1 00 ) 

10 WRITE ( 1 02 v 20) 

20 FORMAT! '.IN THIS BLACKBODY DATA PROGRAM » YOU WILL',/, 
$ 'ENTER A MAGNETUDE AND TEMPERATURE <K ) ' , /, 

$ ' AND THE PROGRAM WILL COMPUTE THE PHOTONS PER',/, 

$ •' WAVELENGTH . " , / , 

$//, 'ENTER THE MAGNETUDE AND TEMPERATURE <K ).' ) 

READ < :l 01 GO.i SM AG, TEMP 
30 FORMAT < 2G ) 


CALCULATE VALUE FOR VISIBLE IRRADI ANCE 

WPC=3. 1E-13*< 10**<-.4*SMAG) ) 

N“0 

VISUAL-0 
DO 50 N==l, 24 

INPUT VALUES OF EYERESPONCE FOR EACH WAVELENGTH 

READ < 20 , 40 ) S < N ) 

40 FORMAT! 1G) 

50 CONTINUE 
N==0 

DO 60 WAVLEN=„ 34 , . S2, . 02 
N--=N-H 

CONVERT MICRONS TO CENTIMETERS 
SD=WAVLEN*lE-4 

CALCULATION OF PLANCK FUNCTION 

W ( N ) =3. 741 85E-12/ ( ( SD*#5 ) * ! EXP ! 1 . 4388/ ! SD*TEMP ) )-l ) ) 


1 
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*TY43~84 

43.000 C 

44.000 C 

45.000 C 

46.000 

47.000 

43.000 

49.000 C 

50.000 C 

51.000 C 

52.000 

53.000 

54.000 

55. 000 

56.000 C 

57.000 C 

53 . 000 C 

59.000 

60. 000 

61.000 C 

62.000 C 

63.000 C 

64.000 C 

65 . 000 

66.000 C 

67.000 C 

63.000 C 

69.000 

70.000 C 

71.000 C 

72.000 C 

73.000 

74.000 C 

75.000 C 

76.000 C 

77.000 

73.000 C 

79.000 

30.000 

31 . 000 

82 . 000 

53. 000 

34.000 


INTEGRATE EYE RESPONSE DATA TIMES PLANCK FUNCTION 

V I SUAL= V I SUAL+W < N ) *S ( N ) *2E~6 
60 CONTINUE 

T0T=5. 679E-12*TEMP**4 

VALUE FOR VISIBLE RESPONCE FRACTION 

EFF=VISUAL/TOT 

N=0 

DO 30 WAVLEN=.34, .32, .02 
N=N+1 


LOCATE PEAK OF BLACKBODY CURVE 

HPEAK= ( WPC/EFF ) * ( ( 1 . 29E- 1 5*TEMP**5 ) / ( 5 . 679E- 1 2*TEMP**4 ) ) 
SA=EXP( 1 . 4330/ < ( WAVLEN#! E— 4 ) #TEMP ) ) -1 

CALCULATE VALUES FOR BLACKBODY CURVE 
AT EACH WAVELENGTH INCREMENT 

H ( N ) = < . 2900/ ( < ( (WAVLEN#1E~4)*TEMP)##5)#SA> ) #HPEAK 
AREA FOR SPACE TELESCOPE EQUALS 3.89E4 CM#*2 
HST (N)=H(N>#3. S9E4* . 02 

CALCULATE NUMBER OF PHOTONS FOR EACH WAVLENGTH 

PHOTON ( N ) =HST ( N ) / ( 1 . 9S64776E-23/ ( WAVLEN* 1 E-4 ) ) 

CONVERT PHOTONS/SEC TO PHOTONS /25M I LL I SEC 

PHOTON ( N ) =PHOTON ( N ) * . 025 

WRITE( 102, 70) WAVLEN, H(N) , PHOTON! N) 

70 FORMAT ( / » 1 X » 'WAVELENGTH" F6 . 3 . 3X , 'H (LAMBDA ) = ',1G, 

*3'X , •- PHOTONS" , F6 . 0 ) 

30 CONTINUE 

STOP' END OF PROGRAM' - 
END 


f 
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Variables in Program and Their Meaning 


SMAG 

Visual magnitude 

Temp 

Temperature of star 

WPG 

Visible irradiance 

S(N) 

EYE response data 

W(N) 

Planck function 

Visual 

Integration value of visible spectrum 

TOT 

Stefan Boltzman function 

EFF 

Blackbody visible response fraction 

HPEAK 

Peak of blackbody curve 

H(N) 

Values for blackbody curve 

Hst(N) 

Values represented in watts 


Photon (N) 


Number of photons at each wavelength 



LIST OF EYE RESPONSE DATA FOR PROGRAM 


! E EYERESP 
EDIT HERE 
*TY 


.340 0.0 
.360 0.0 
.380 0.00004 
.400 0.0004 
.420 0.0040 
.440 0.023 
.460 0 . 060 / 
.480 0. 139 
.500 0.323 
.520 0.710 
.540 0.954 
.560 0.995 
.580 0.870 
.600 0.631 
.620 0.381 
.640 0. 175 
.660 0.061 
.680 0.017 
.700 0.0041 
.720 0.0 
.740 0.0 
.760 0.0 
.780 0.0 
.800 0.0 
.820 0.0 

* 


f 
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